Memory-keeping effects and forgetfulness in the dynamics of a qubit coupled to a spin chain 
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Using recently proposed measures for non-Markovianity [H. P. Breuer, E. M. Laine, and J. Piilo, Phys. Rev. 
Lett. 103, 210401 (2009)], we study the dynamics of a qubit coupled to a spin environment via an energy- 
exchange mechanism. We show the existence of a point, in the parameter space of the system, where the 
qubit dynamics is effectively Markovian and that such a point separates two regions with completely different 
dynamical behaviors. Indeed, our study demonstrates that the qubit evolution can in principle be tuned from a 
perfectly forgetful one to a deep non-Markovian regime where the qubit is strongly affected by the dynamical 
back-action of the environmental spins. By means of a theoretical quantum process tomography analysis, we 
provide a complete and intuitive characterization of the qubit channel. 

PACS numbers: 03.65. Yz, 75.10.Pq , 42.50.Lc 



43 



> 
CO 

m 



In the study of open quantum systems, quite often the 
Markovian approximation has been a useful starting point to 
describe their dynamics. Only very recently, new and power- 
ful tools have been designed in order to tackle the important 
question of explicitly quantifying the non-Markovian charac- 
ter of a system-environment interaction or a dynamical map 
Ql-01 • Such a task is extremely important, given that non- 
Markovian effects are known or expected to occur in a wide 
range of physical situations, especially in the realm of solid- 
state devices where a system of interest is often exposed to 
memory-preserving environmental mechanisms. As a specific 
example, switching impurities have been shown to affect su- 
perconducting devices in various regimes 0]. 

From a fundamental point of view, devising a reliable way 
to actually quantify non-Markovianity is very useful in light 
of the plethora of frequently ad hoc or technically rather 
involved approaches put forward so far in order to study 
memory-keeping environmental actions. A few instances 
emerge as promising measures for the non-Markovian nature. 
Wolf et al. IU\ have proposed to quantify the degree of non- 
Markovianity of a map by considering the minimum amount 
of noise required in order to make the evolution of a sys- 
tem fully Markovian. On the other hand, in Ref. |2], Rivas 
and co-investigators discussed two approaches, founded on 
the deviations of a given map from full divisibility. Finally, 
Breuer et al. proposed a way to quantify non-Markovian ef- 
fects by looking at the back-action induced on the system 
under scrutiny by its memory-keeping environment yfl. Fur- 
ther approaches have been considered and measures have been 
proposed which are based on the use of other interesting in- 
struments such as the quantum Fisher information (5D. 

In this paper, we focus on the non-Markovianity measure 
proposed in Ref. yD to study the dynamics of a qubit coupled 
to a spin environment described by an XY model in a trans- 
verse magnetic field. Our aim is to analyze a simple and yet 
non-trivial system-environment set, displaying a broad range 
of behaviors in the parameter space, in order to relate the fea- 
tures of the non-Markovianity measure to the known spectral 
and dynamical structures of the spin system. The choice of 



the model has been also dictated by the versatility demon- 
strated by spin-network systems in the engineering of proto- 
cols for short-range communication [6] and also for the in- 
vestigation of the interplay between quantum- statistical and 
quantum-information-related aspects ED. 

Intuitively, one would expect that the coupling to a spin 
environment always leads to a non-Markovian dynamics for 
the qubit. This can be justified by noticing that the environ- 
mental correlation time is non-zero even in the weak coupling 
regime and by conjecturing that, by increasing the coupling 
strength, a Markovian approximation can only become less 
valid. However, we demonstrate that this is not at all the 
case and that peculiar behaviors occur at intermediate cou- 
plings. The chosen measure for non-Markovianity, in fact, 
turns out to be identically null at specific dynamical regimes, 
thus demonstrating the absence of a net re-flux of information 
from the environment back to the qubit. Our study of non- 
Markovianity is thus novel and makes the use of measures 
for its quantification a valuable tool for the identification of 
special conditions in the parameter space of the environment. 
Moreover, our study opens up the possibility to exploit the 
different memory-keeping regimes induced by controlling and 
tuning the properties of the environment to effectively "drive" 
the qubit dynamics in a non-trivial and potentially very inter- 
esting way. The effective evolution of our two-level system 
could be guided across different regimes, ranging from strong 
environmental back-action to completely forgetful dynamics 
typical of a Markovian map. The potential of such flexibility 
for reliable control at the quantum level is a topic that will be 
explored in future. 

The reminder of this paper is organized as follows. In Sec.U 
we introduce the model under scrutiny and revise the basic 
principles behind the chosen measure for non-Markovianity. 
We study a few cases amenable to a full analytical solution 
and highlight how, in a few of such instances, a simple ex- 
perimental protocol can be designed for the inference of the 
properties of the qubit dynamics. Most importantly, we re- 
veal the existence of an operating regime where the chosen 
measure is strictly null and relate such an effect to intriguing 
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modifications occurring at the level of the energy spectrum of 
the qubit-environment system. Sec. HU is devoted to the for- 
mal characterization of such peculiar point in the parameter 
space. We first demonstrate divisibility of the corresponding 
dynamical map, hence its Markovian nature, and then perform 
a theoretical analysis based on the use of quantum process 
tomography to quantitatively infer its properties. Finally, in 
Sec. ITU we draw our conclusions. 



I. THE MODEL AND THE MEASURE 

We consider a qubit Q coupled to a chain r of N interact- 
ing spin- 1/2 particles. The qubit is described by the spin- 1/2 
vector operator s > while the operator s n (n=l, .., N) corre- 
sponds to the spin located at site n of the chain Y. The log- 
ical basis for the spins and Q is given by {|0)j , |l)j} with 
j = 0, .., N. The Hamiltonian ruling the intra-chain interac- 
tion is taken to be of the XF-Heisenberg type (we set H = 1) 
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where h n is the local field applied at site n. T is open-ended 
with J x,y 's and h n 's being not necessarily uniform. The qubit 
is coupled to the first spin of the environment, embodied by T, 
via an exchange interaction of strengths Jq ,y and is subjected 
to a local field ho according to 



7} = -2(«^ + J y sls\)-2h sl. 



(2) 



In order to determine the time-evolution of Q, we resort to 
the Heisenberg picture and the formal apparatus put forward 
in Refs. Ill Oil , which provides particularly powerful tools for 
the study of the many-body problem embodied by Ho + l-Lv • 
Using the operator-expansion theorem and the algebra satis- 
fied by Pauli matrices, we find that the time evolution of the 
components of so reads 



(3) 



n=0 
1 N 

sl{t) = -Y J [Kl{t)vi-K(t)K\Pn, 

s z o(t) = -ia* (t)&Ut)/2, 



where a% (a=x, y, z) are the Pauli operators for the spin 
at site n and P n = Yl^i °t • The time-dependent coeffi- 
cients and A x (t) are the components of the (N + 1)- 
dimensional vectors U x (t) and A x (t) defined by 

n*(t) = E(- 1 ) p ttit( ttT )^' ( 4 ) 

= E(-!) P ^ I ^, t T (tt t Tv, (5) 



where T stands for transposition, the vector v has components 
V{ = Sio and the tri-diagonal adjacency matrix r has elements 
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Notice that we have labelled columns and rows of 
(iV+l)x(iV+l) matrices and (iV+1) -dimensional vectors 
using indices ranging from to TV. The coefficients Yi y (t) 
and A y (t) are obtained from Eqs. © and (0) by replacing r 
with r T . Both rr T and r T r can be easily diagonalized by 
orthogonal matrices U and V such that (rr T ) p =UA 2p U T 
and (r T r) p =VA 2p V T , with A a diagonal matrix whose ele- 
ments Xij = XiSij are the (positive) square roots of the eigen- 
values of t t t. Consequently, Eqs. (01) and (0) can be fully 
summed up to give 



U x (t) = Un(t)U T v, A x (t) = W(t)U T v, 
U y (t) = VQ(t)V T v, A y (t) = UV(t)V T v, 



(7) 



where fl(t) and are diagonal matrices with elements 
Qij(t)= cos(\it)5ij and E^-(£)= sin(A^)^j. 

By using Eqs. ([3]), one can determine the time evolution of 
the state of Q as p(t)=i/2+ Ja^o C0)^o • In the evalua- 
tion of the expectation values required to determine p(t) we 
assume that Q and V are initially uncorrelated. The conserva- 
tion rule [H, ®n=o °n\ = and me property (^n^m^n) =0 
imply 



n=0 
1 N 

W)>=2 E [WK^> - A y n (t)(P n a x ) 

1 N 

<*o(*))=2 E KW^W + K(t)A y n (t)} (K) 

n=0 



p=0 



(2p+l)! 



n<m 

-- $3 K(t)n^(t) + A^(t)A^(t)] (p„ +1 p m ^>. 

n < rn 

(8) 

In order to evaluate the above equations one needs multi- 
spin correlation functions, involving, in particular, the degrees 
of freedom of T. To this end, we consider it to be in its 
ground state in what follows. However, as it will be shown 
later on, the value of the measure of non-Markovianity cho- 
sen for this work is independent of the state of T, provided 

that (s?^ (0)) =0 holds, as it does for the ground state with 
no broken symmetry. 

By using Eqs. ([5]), one can finally determine p(t). The in- 
teraction with the spin chain acts for Q as a dynamical map 
$(t,0) such that p(t) = $(t,0)p(0). The properties of the 
map depend on the relative weight of the various parame- 
ters entering tio+tir- Our aim is to characterize the non- 
Markovian nature of as a function of these parameters. 
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To pursue our task, we consider the measure proposed in 
Ref. OD, based on the study of the time-behavior of the trace 
distance V[p^\t), p (2) W]=|Tr|pW(t)-p( 2 )(t)|G[0, 1] be- 
tween two single-qubit density matrices p^ 1,2 \t). The trace 
distance is such that V[p^\t) ^ p^ 2 \t)]=l when the two 
probed states are completely distinguishable, while it gives 
for identical states 111311 . The degree of non-Markovianity 
Af(<&) of the dynamical map <£, is defined as 

JV($)= max^{X>[pW (b n ), (b n )] -X>[p« (a n ), p^ (a n )} } 

n 

(9) 

where the maximization is performed over the states p^ 1,2 ^ (0) 
and (a n , b n ) is the n th time window such that 

a[t,pW(0)] = d t V[pV(t),pM(t)] > 0. (10) 

The function cr[t, p( 1>2 ) (0)], which has been dubbed^wx of in- 
formation in Ref. [3], encompasses per se the condition for 
revealing non-Markovianity of an evolution: the mere exis- 
tence of even a single region where cr[t, p( 1,2 ) (0)]>0 is suffi- 
cient to guarantee the non-Markovian nature. Conceptually, 
in fact, Af($) accounts for all the temporal regions where 
the distance between two arbitrary input states increases, thus 
witnessing a re-flux of information from the environment to 
the system under scrutiny. Such re-flux of information ampli- 
fies the difference between two arbitrarily picked input states 
evolved up to the same instant of time. A Markovian dy- 
namics is such that the above-mentioned re-flux never oc- 
curs and a[t, p^ 1,2 \0)} < always. For the case at hand 
and for two generic input density matrices p( 1,2 \0), we find 
D[pW(t),pM(t)} = VC where 



and the environmental chain T as 
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C=PS(t)IIg(t) + AS(t)Ag(t)]V 
+ | c/+ (IIo , Ao , t) + c* /_ (IIo , A , t) 1 2 /4 



(11) 



with/±(n , A ,t)=IIg(t) ±ng(t)±2 [Ag(t)±Ag(t)] where 
P=PnC)-Pn(0) and c=p^ ) (0)-p^ ) (0). We have used 
the notation p[j\o)= (i\p^ \ j) Q with i,j=0, 1 and a=l, 2. 
It is worth noticing that in the above equation the initial state 
of the environment is completely absent, so that the environ- 
ment's multi-spin correlators are not relevant. Eq. (TTTh can 
be recast into a much more intuitive form by referring to the 
Bloch vectors r a =(r^, r^, r z a ) representative of the state p^ 
(we have chosen r a = (sq) as the mean value of the spin- 
1/2 operator rather than that of the more usual Pauli operator 
to avoid the appearance of irrelevant 1/2 factors). By calling 
^ r a =r a_ r a difference between the vectors of two input 
states, we have 

C=|r 1 (t)-r 2 (t)| 2 ^[ng(t)ng(t)+A§(t)Ag(t)] 2 (Ar^) 2 

+ [n x (t)Ar x +A x (t)Ar y } 2 + [U y (t)Ar y -A v (t)Ar x ] 2 . 

(12) 

One can thus write the flux of information between the qubit 



"[*,P (1 ' 2) (0)] 



(Ar z Y Ad t A + BdtB + Cd t C 



(13) 



where A = U x {t)U y (t) + Ag(t)Ag(t), B = U^(t)Ar x + 
Ag(t)Ar^, and C = U y (t)Ar y - A y (t)Ar x . Due 
to the non-negativity of V[p^(t), p^(t)], the condition 
for non-Markovian dynamics can be simply stated as 
(Ar z ) 2 Ad t A + Bd t B + Cd t C>0. Numerically, it turns out 
that the maximum in the corresponding measure of non- 
Markovianity is achieved for p( 1,2 \0) being antipodal pure 
states lying on the equatorial plane of Q's Bloch sphere (we 
will come back to this point later in this paper). This consid- 
erably simplifies the necessary condition for memory-keeping 
dynamics to the form d t (B 2 +C 2 )>0. 

Although our analysis can be carried out without major 
complications in the general case, in order to simplify the pre- 
sentation, from now on we restrict ourselves to the case of a 
uniform spin environment with equal isotropic couplings be- 
tween every pairs of nearest neighboring spins (XX -model) 
and set this coupling constant as our energy (and inverse time) 
unit. We thus consider Jq = Jq= Jo, a condition under which 
the measure of non-Markovianity becomes 



V[pM(t),p™(t)] = J(p*f(t) + \c\ 2 )f(t), 



(14) 



where f(t) = II§ (t) + A§(t) € [0, 1]. The corresponding rate 
of change of the trace distance is 



*[*,P (1 ' 2) (0)] 



(2p 2 /(t) 



fit) 



2V[p^)(t\p^{t)] 



(15) 



As f(i)>0, the sign of a[t, p (1 ' 2) (0)] is determined by f'(t), 
regardless of the pair of input density matrices p^ 1 ' 2 ) (0). Yet, 
we should look for the states maximizing the contributions to 
the trace distance within these time intervals. As |/(t)|<l, 
such optimization is achieved for c=l and p=0; that is, an- 
tipodal pure states on the equatorial plane of the Bloch sphere, 
in line with the (more general) numerical findings reported 
above. The condition for non-Markovianity can be further 
elaborated as f (t)=2[U (t)U f (t) + A (t)A 7 (£)]>0. 

To pursue the task of evaluating A/"($) in the quite rich 
parameter space of our model, we start by considering the 
case of a qubit resonant with the spin environment, and as- 
sess the special case of ho=h first. It can be shown analyt- 
ically that in this case one gets Unit) = J\{2t) cos(2ht) / 1 
and A (t) = -Ji(2t) s'm(2ht)/t |9], where J n (x) is the 
Bessel function of order n and argument x. Correspondingly, 
f(t)=J 2 (2t)/t 2 and/ 7 (t) = -4Ji(2t) J 2 (2t)/t 2 . For those 
states with |c| = 1 and p = which maximize J\f, we have 
the flux 

a[t,p^ 3 \0)} = -(2/t)spx[J 1 (2t)]J 2 (2t) (16) 

with sgn[x] being the sign function. Eq. (TT61) is independent of 
h, which is joint result of the condition p = and the invari- 
ance of the trace distance under the global unitary transfer- 
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FIG. 1: (Color online) A/"($) versus h/J with h /J = 0,J /J=1 
and N — 100. To avoid spurious recursion effects, the dynamics is 
evaluated up to a temporal cut-off of ~ 2N/3. We checked that the 
precise cut-off is not relevant and that the plot remains unchanged if 
N is varied, provided N ^> 1. All quantities are dimensionless. 



mation embodied by the operator exp[— ith Y^iLo °V\ nee ded 
in order to pass to the interaction picture. The time windows 
where a[t, p^ 1,2 ^ (0)]>0 are determined by the chain-rule of 
the zeros of the Bessel functions. Overall, we get that the flux 
is positive for a time t a> o=U(tj,t 2 ) where t\ (t 2 ) is the i th 
zero of the Bessel function of order 1 (2). From this special 
case and the general considerations reported above, we learn 
that A/^) only depends on the detuning Sh = h — ho. 

Due to this fact, the simple case in which all of the mag- 
netic fields are absent (h = ho = 0), which can be tested 
experimentally in an easy way, allows to draw some interest- 
ing conclusions on the more general case Sh = 0. In absence 
of magnetic fields, indeed, the behavior of a and M can be 
extracted directly by monitoring the dynamics of the qubit. 
To show that this is indeed the case we start by noticing that, 
if h/J = h 1 J = 0, we have A (t)=0 |9], so that both 
the trace distance and the flux of information are determined 
by f(t) = IIq (t), which represents the squared length of the 
Bloch vector, \r(t)\ 2 . This implies that, in order for Q to ex- 
perience a re-flux of information from T, IIq(£) has to be a 
strictly non-monotonic function. Therefore, we observe non- 
Markovian dynamics when the Bloch vector of the qubit is 
alternatively shrunk and elongated during its evolution. This 
can be witnessed by reconstructing the density matrix of Q 
using standard quantum state tomography techniques fj"ll [l3ll. 
which are routinely implemented in a variety of physical se- 
tups. However, there is also an interesting alternative which 
does not require full state reconstruction. From the first of 
Eqs. ©, we get that (sg (t)) = §n (t)(£jj (0Y). Thus, in the 
spirit of the proposals put forward in Refs. Oil [12], by prepar- 
ing the state of Q in one of the eigenstates of <Tq [a choice 
that would be perfectly consistent with our results on the input 
states to be used for the calculation of A/"(<£)] we can measure 
(so(t)) to determine the non-Markovianity of the qubit dy- 
namics, which would be revealed by its non-monotonic time 
behavior. 

Other points in parameter space exist for which the model 
is amenable to an exact analytic solution. However, in order 
to give a complete overview of the behavior of Af, we first 



resort to numerical techniques to solve our model. The results 
of such an analysis are shown in Fig.Q] where the quantitative 
degree of non-Markovianity is shown as a function of h/J 
(we set ho=0) in the isotropic case and for equal intra-chain 
and qubit-environment coupling strengths. 

A highly non-trivial behavior followed by A/"(<£) is re- 
vealed. The largest deviation from a Markovian dynamics is 
achieved at h/J = while, for < h/J < 0.5, A/"($) de- 
creases monotonically to zero. By further increasing h we see 
that Af(&) achieves a very broad maximum around the sat- 
uration point h/J ~ 1 of the environmental chain. Finally, 
it goes to zero for h/ J ^> 1, as it should be expected given 
that this situation corresponds to an effective decoupling of the 
qubit from the environment [9] . By generalizing our study to 
the case of ho/ J ^ 0, we emphasize the presence of a Marko- 
vianity point at Sh/ J =1/2. It turns out that this point sepa- 
rates two regions in which the dynamics of the qubit (although 
being non-Markovian in both cases) is completely different. 
For Sh/J<l/2, indeed, the qubit tends toward a unique equi- 
librium state at long times, irrespective of the initial condition; 
that is, the trace distance goes to zero after some oscillations. 
For larger detunings, on the other hand, the trace distance does 
not decay to zero, implying that some information about the 
initial state (and in particular, about the relative phase between 
its two components) is trapped in the qubit. 

The Markovianity point at Sh/ J =1/2 (with 
Jo/J=l) is one of those points in parameter space 
for which the model can be treated fully analyti- 
cally Q. We have U (t)=J (2t) cost+Ji(2t) sint 
and A (t)=Ji(2i) cost- J (2t) suit Jl, so 
f(t)=J 2 (2t)+J 2 (2t) and f(t)=-2J 2 (2t)/t<0 V*. 
As a consequence, the corresponding non-Markovianity 
measure is always zero at this point of the parameter space, 
which is a very interesting result due to the relatively large 
Q — T coupling strength. Such a feature, in fact, would 
intuitively lead to exclude any possibility of a forgetful 
dynamics undergone by the qubit. Yet, this is not the case and 
a fully Markovian evolution is in order under these working 
conditions. A deeper characterization of this Markovian 
dynamical map is given in Sec.JIl 
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FIG. 2: (Color online) JV($) versus h/J for h o /J = and three 
different values of Jo/ J for a chain of N = 200 spins. For Jo / J> 1, 
the Markovianity point disappears. All quantities are dimensionless. 
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As Jo represents the energy scale of the qubit-environment 
interaction, it is natural to expect that significant changes in 
Af(&) occur as this parameter varies. In particular, we find 
that i) a Markovianity point with J\f(<f>)=0 only exists for 
Jo/J<l and ii) Af($) tends progressively towards a mono- 
tonically decreasing function of Sh if Jo/ J grows from 1 
to y/2 [see Fig. 0. Strikingly, at J /J=v / 2 and 5h=0, 
the adopted measure of non-Markovianity diverges, as it can 
be checked by using the analytic integrability of the qubit- 
chain interaction at this point in the parameter space. In- 
deed, at J /J=v / 2, we have U (t) = J (2t) and cr(t) = - 
2sgn[j7o(2£)]^7i(2£). Integrating over all the positive time in- 
tervals determined by means of the usual chain rule we get 
A/"($)— )>oo, thus witnessing a strong back-action of T on the 
state of the qubit. The divergence of Af(&) should not surprise 
as it is common to other situations with spin-environments, 
such as the so-called central- spin model where a single qubit 
is simultaneously coupled to N independent environmental 
spins via Ising-like interactions (see Breuer et al. in Ref. yfl). 
We provide a physical explanation for the enhanced non- 
Markovian nature of the qubit dynamics simply by looking 
at the spectrum of the Hamiltonian ruling the evolution of the 
qubit-chain system. For Jo / J>V%, the spectrum of tir + %o 
exhibits a continuous spectrum (a band of extended eigen- 
states) that is lower- and upper-bounded by two discrete en- 
ergy levels whose eigenstates are localized at the sites occu- 
pied by the qubit and the first spin of T [15]. As a conse- 
quence, a certain amount of information remains trapped into 
such a localized state, bouncing back and forth between the 
qubit and the first spin and therefore mimicking a highly non- 
Markovian dynamics characterized by strong back-action, so 
that Af($) diverges. 

An analysis similar to the one performed just above allows 
us to obtain an intuition for the behavior of the measure of 
non-Markovianity near this point. For Sh/J>l/2 the spec- 
trum of the system shows one eigenenergy out of the band 
and its corresponding eigenvector is localized around the site 
occupied by Q. This can explain the information trapping that 
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FIG. 3: (Color online) We call E+ (E-) sl discrete level lying above 
(below) the continuous energy band in the (single particle) spec- 
trum of the total Hamiltonian Tio + Hr- The parabolae h/J = 
±[1 — (Jo/ J) 2 / 2] divide the (h/J, Jo/ J) parameter plane in three 
regions with one, two and no localized state (light-colored, dark- 
colored and uncolored region in the plane, respectively). All quanti- 
ties are dimensionless. 



occurs for Sh/ J>l/2. For definiteness, in what follows we 
report explicit results for the case ho/ J = 0. Therefore, from 
now on, we consider 5h/J = h/J. 

In order to determine the existence of a more general con- 
nection between the emergence of localized eigenstates in the 
system-environment spectrum and a point of zero-A/" in the 
qubit evolution, we analyze the (h, Jo) -plane to find out where 
localized eigenstates appear and then evaluate the correspond- 
ing degree of non-Markovianity. In doing this, we take ad- 
vantage of the fact that the Hamiltonian describing an en- 
vironmental XX-model has the same single-particle energy 
spectrum as a tight-binding model with an impurity. Follow- 
ing the approach given in Ref. |16], we deduce that, in the 
(h/J, J / J) plane, the parabolae h/J = ±[l - (J /J) 2 / 2 ] 
define regions with respectively zero, one and two localized 
energy levels out of a continuous-energy band (see Fig.0. 

The central region with no localized energy state and the 
zones with only one localized state (either a upper-lying or 
lower-lying one with respect to the continuous band) cor- 
respond to Af(<f>) ^ 0, although finite. The frontiers of 
such regions, marked by the parabolae, give the limiting val- 
ues of (h/J, Jo/ J) for the appearance of the localized states 
and are such that Af(&) = only for Jo /J < 1, whereas 
for 1 < Jo / J < y/2 the measure of non-Markovianity is 
non-null and finite. Finally, the regions with two localized 
states and their frontier with the no-localized state region have 
Af($)—>oo, in line with the analytical results discussed above. 

In particular, we remark once again that the Markovian- 
ity points (when they exist, i.e. for J / J < 1) stay on the 
parabolae; that is, they are found to occur at the onset for the 
existence of one discrete eigenstate outside the energy band. 
This kind of state contains a spatially localized spin excita- 
tion with a localization length that decreases with increasing h 
Jl5ll . Precisely at the border, the localization length becomes 
as large as the length of the chain itself, so that all of the en- 
vironmental spins are involved in (or share excitation of) the 
initial state. This can intuitively justify the fact that J\f is zero 
there. For Jo /J > 1, on the other hand, the first spin of the 
environment becomes more important for the dynamics of the 
qubit (which is "strongly coupled" to it) and an information 
exchange is always found to occur between them, irrespec- 
tively of the existence of the discrete level. This information 
re-flux becomes more and more pronounced with increasing 
Jo / J and decreasing h/J. 

To stress once more the close relationship between the non- 
Markovianity measure and the properties of the overall Hamil- 
tonian 1-Lt we report here the energy distribution of the 
excitations which are present in the initial state of the system 
(given by the product of an equatorial state for the qubit times 
the ground state of the environment). These excitations are 
spin-less fermions of the Jordan- Wigner type and the proce- 
dure to obtain them is the one described, e.g., in Ref. lfl7h . 
Fig. [4] reports the average value of the excitation number 
in the initial state vs single-particle-energy for various val- 
ues of the magnetic field. The plot shows what happens near 
a Markovianity point: the energy distribution of initial-state 
excitations becomes flat (i.e., structure-less) at the J\f = 
point, while it shows a maximum for lower values of h, cor- 
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responding to a finite value of J\f, and a spike for h/ J > 0.5, 
corresponding to the discrete level giving rise to information 
trapping. 

From this discussion we conclude that the measure of non- 
Markovianity of the qubit dynamics is in fact a detector of 
general aspects of the full qubit+environment system and that 
various features of J\f(<&) can be related to general character- 
istics of the overall many-body problem described by the full 
Hamiltonian model. 



II. CHARACTERIZATION OF THE POINT OF 
ZERO-MEASURE 

The occurrence of a null value of Af(&) ath/J= 1/2 and 
Jo/ J — 1 deserves a special attention. Naively, one could ex- 
pect a special behavior to occur at the chain saturation point 
(i.e. at h/J = 1), where the intrinsic properties of the envi- 
ronmental system are markedly different from the situation at 
h/J < 1. To the best of our knowledge, indeed, no signifi- 
cant dynamical feature has been reported for the model under 
scrutiny away from saturation. 



A. Characterization of the dynamical map: formal features 
and divisibility 

The aim of this Section is to characterize the dynamical 
map that we obtained for the qubit under these conditions. 
Thanks to the analytical solution for the case at hand provided 
in Ref. fl4ll . we can sum up the terms appearing in Eq. (8]) and 
determine the complete density matrix of Q. For n = 0, .., N 
we have 



n n (t)= 



(-l)5[sin(t) t 7 n+ i(2t)+cos(t)J r n (2t)] 
(-1)" ' 



n even, 



[C0s(t)j r n +i(2t)- 



sm(t)J n (2t)] nodd 
(17) 



0.06 



0.04 



0.02 



n(k) 




energy 



1 



FIG. 4: (Color online) Average number of fermion excitations in the 
initial state for a qubit-chain system with N + 1 = 50, coupling 
strength Jo/ J = 1 and three values of magnetic field around the 
Markovianity point (h/J — 0.5, ho /J = 0). The initial state is 
always taken to be the tensor product of an equatorial state of Q and 
the ground state of T. 



and 



A (t)= J (-l) ? [cos(t)J n+ i(2t)-sin(t)J n (2t)] neven, 
nU \(-l) I ^[sin(t)J n+1 (2t)+cos(t)J n (2t)] nodd. 

(18) 

The dynamical map transforms the elements of the input den- 
sity matrix pij(0) (i, j = 0, 1) as 

Poo(t) = /(t)poo(O) + Aj (t), p 01 (t) = Agl(t)poi(0), 

pii(t) = + /(f)pn(0), P io(t) = iOt)pin(O), 

(19) 

= [l-f(t)-g(t))/2 
and where 
defined in Sec. [H 

= ELiKW + 

A n (t)A m (t)]g 

nm-> 

and the two-point 



with AH{t)=[l-f{t)+g{t)]/2, A\\{t) 
and A° \(t)=Ilo(t)+iA (t)=[A{°(t)Y 
f(t)=U (t) 2 +A (t) 2 has already been 
We introduced here the function g(t) 

A^(t)](^>-jEL m =i[nn(t)n ro (t)4 

written in terms of the magnetization 
longitudinal correlation function 



<*n> = l- 



N+ 1 

Qnrn = (J^n Prn & n & rn ) = 



cos[(/cf + l)l? n ] sm[kF$ r , 
sin $ n 



2 (cOS$ n — COS^ m ) 



(20) 



with ^■ >ib = v /2/(JV + l)sin(ji? fe ), = ^/(7V + 1), 

fc = l,..,iV and fci? being the Fermi wave number, see 
Ref. 



The condition for divisibility stated in Ref. [3] implies the 
existence of a completely positive dynamical map ^f(t + 1\ , t) 
such that, for two arbitrary instants of time t and ti, we have 
$(t + ti,0) = *(t + ti,t)$(t,0). Here$(t,0) is the dynam- 
ical map in Eq. (fT9l ). Any dynamical map that is divisible ac- 
cording to the above definition is Markovian. This implies that 
non-divisibility is a necessary condition for memory-keeping 
effects in the evolution of a system. A dynamical connection 
t) between the states 0)p(0) and 0)p(0) 

can be straightforwardly found to be given by the map chang- 




FIG. 5: (Color online) Divisibility condition C against t and t\ for 
qubit Q initially prepared in (|0) + |l) )/v / 2- Divisibility is guar- 
anteed for < C < 1. In panel (a) we have taken Jo/ J = 1 with 
ho/ J = h/J = 1/2, while in (b) it is h/ J = 1.1. All quantities are 
dimensionless. 
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ing the elements of the qubit state p(t) at time t into 



with it. 



/(* + *i) 
/(*) 



oo„ , ^ ^ + ^00^ 



M* + 0=P*o(* + *i) 



/(*) 



(21) 

Therefore, in order to ensure the divisibility of <!>(£ + 
ti,0), we should investigate the complete positivity of 
+ ti,t). To this purpose we make use of the Choi- 
Jamiolkowski isomorphism 111 811 and prove the complete pos- 
itivity of + t±,i) by checking the non-negativity of 
(%f(t + t\, t) (8) I2) P, where p is the density matrix of one of 
the Bell states |19] and I2 is the identity map. By choos- 
ing p(t)=\ (|00) + |11)) (<00| + (11|), the action of the map 
(%f(t + ti, t) I2) determines the following non-zero matrix 
elements (up to an irrelevant factor 1/2) at time t + ti 

Poo,oo(t + fx, t)=^±M+AjJ(t + tO-^y^^SW, 

Poo,n(* + ti,t)=AH(t + tl )- f ^ t ^ A&(t), 

00 u , ^ / ^ + tl ^??(t), 



Pn,oo(* + *i,*)=AS(* + *i) 



/(*) 



Pll,ll(^ + = 



/(*) 



+ASS(t + ti)-^^A?S(t), 



/(*) 

' 



A)i,oi(* + *i,*)=Poi 01 + *1j *) 

(22) 

Here, we have used Pi>j>,ij= \ p with the primed (un- 
primed) indeces corresponding to the evolving (non-evolving) 
qubit. The condition for positivity of the composite two-qubit 
map turns out to be equivalent to the positivity condition of the 
single-qubit one given in Eqs. (I2TK which is in turn translated 
into the inequality < C{t,t\) < 1 with 

C(Mi) = 4|poi(^ + ^i)| 2 + [poo(t + ti)-pn(t + ti)] 2 . (23) 

In Fig. [5] we show the typical behavior of C(Mi) at me 
Markovianity point [panel (a)] and away from it [panel (b)] . 
We have taken the qubit Q as prepared in (|0) + \l) )/V2, 
which is a significant case as equatorial states in the Bloch 
sphere are those optimizing the calculation of J\f(®). Al- 
though this is simply a representative case, we have checked 
that for a uniform distribution of random initial states of Q, 
no significant quantitative deviations from the picture drawn 
here are observed. Clearly, by moving away from h/J=l/2, 
temporal regions where C(£, ti)>l are achieved. This demon- 
strates, from a slightly different perspective, the flexibility 
of the effective qubit evolution: a wide range of dynami- 
cal situations is spanned, from fully forgetful to deeply non- 
Markovian dynamics, strongly affected by the environmental 
back-action. The kind of evolution of Q can be determined by 
tuning the parameters of the environment and its interaction 



B. Formal characterization of the channel through theoretical 
quantum process tomography 

We now turn our attention towards the formal characteri- 
zation of the channel achieved at the Markovianity point, so 
as to qualitatively explain the reasons behind the nature of 
the corresponding qubit evolution. We stress that this sort of 
investigation is meaningful only at this specific point in pa- 
rameter space, where complete positivity is guaranteed. In 
principle, full information on the reduced dynamics of Q 
could be gathered from the Kraus operators {Ki} such that 
0)p= J2i KiqK\ w i m P m e density matrix of the generic 
qubit and KjKi=i. As their direct computation is not 
possible due to the complications of the Q — T coupling, here 
we gain useful information on the structure of the i\Ys by 
means of the formal apparatus for quantum process tomogra- 
phy Il3ll . that we briefly remind here. 

The characterization of a dynamical map reduces to 
the determination of a complete set of orthogonal opera- 
tors {JC m } over which one can perform the decomposition 

K t = e imKm SO as to get 



(24) 



where the channel matrix Xmn= e im^* n has been in- 
troduced. This is a pragmatically very useful result as it 
shows that it is sufficient to consider a fixed set of opera- 
tors, whose knowledge is enough to characterize a channel 
through the matrix \- The action of over a generic ele- 
ment \n)(m\ of a basis in the space of the 2 x 2 matrices (with 
n, m=0, 1) can be determined from a knowledge of the map 
$ on the fixed set of states |0) , |1) , |+) =(1/a?2)(|0) + |1)) 
and 1+3,) =(l/ v / 2)(|0) +i |1)) as follows 

£(\n) (m\)=£{\+)(+\)+i£{\+y)(+y\) 

-^[£(\n)(n\)+£(\m)(m\)]. ^ 

Therefore, each qj = \n) (m\ (with j = 1, .., 4) can be found 
completely via state tomography of just four fixed states. 
Clearly, &(Qj) = J2k ^jkQk as {gk} form a basis. From the 
above discussion we have 

m,n m,n,k k 

where we have defined JCmQjJC^ = ^ k P^Qk so that we 
can write 



(27) 



The complex tensor /3J^ n is set once we make a choice for 
{JCi} and the Aj^'s are determined from a knowledge of 
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®(Qj). By inverting Eq. d27lh we can determine the channel 
matrix x and characterize the map. Let W be the operator di- 
agonalizing the channel matrix. It is straightforward to prove 
that, if Di are the elements of the diagonal matrix W'xjU, then 
e irn = ^/DlUmi so that 

Ki = y/Dl^UjiKlj. (28) 

3 

This apparatus can be applied to the problem at hand. To 
this end, we numerically determine the evolution of the four 
qubit states |0), |1), |+) and \+ y ) at ft/J=l/2, h =0 and 
in the uniform-coupling case for N +1=150. As the time- 
scale within which an excitation travels back and forth to Q 
scales as N BSD, in order to avoid recurrence effects we limit 
the time- window of our analysis to [0, 27V/3]=[0, 100], which 
we chop into small intervals of amplitude 0.2. At each instant 
of time in such a temporal partition, we evaluate the evolution 
of the four probing states given above. This is the basis for 
the reconstruction of the Kraus operators, which is performed 
as described above. The results of our numerical calculations 
are four Kraus operators, whose form depends on the instant 
of time at which the dynamics is evaluated. In principle, an 
analytic form of such operators can be given. Upon inspec- 
tion, one can see that the action of the channel embodied by T 
over the four probing qubit states can be summarized as fol- 
lows 

,„v fa \ | , v fa b + ic\ 

\0)^eo=[p i_ a J.I+>-*+=( 6 _fc i_ J. 

ii \ f P \ | , , / a —c — ib\ 

^^={0 l-p)i\+y)^y={- c + ih 1-a) 

(29) 

with a, a, f3 G [0, 1] and 6, c G C. The form of the cor- 
responding analytical Kraus operator can then be found, al- 
though their cumbersome nature makes them unsuitable to be 
reported here. Despite the lack of analytic formulae, useful 
information can be gathered from this study. Fig. [6] shows the 
dynamics of 10 random pure input states of the qubit. Each 
dot represents the state of the qubit at time t in the space of the 
density matrices. Colors identify the evolved states of a given 
input one. As time increases, i.e. as we follow the points of 
the same color from the outer to the inner region of the plot, 
the Bloch vector of each initial pure state spirals converging 
towards the center. That is, the state becomes close to a max- 
imally mixed state. 

We should mention that, for any fixed N, this does not ex- 
actly occur as the final state of Q lies slightly off- set along 
the z-axis and does not describe a perfect statistical mixture. 
We ascribe such an off- set to the finite length of the chain 
being studied and the existence of a sort of effective magnetic 
field (scaling as 1/N) acting on Q and induced by its coupling 
with the environmental chain, which determines the residual 
polarization of the qubit. Numerical evidences suggest that, 
in the thermodynamical limit, such a residual polarization ef- 
fectively disappears and the final state of the evolution is 1/2. 

As discussed in the previous Section, from the study of 




FIG. 6: (Color online) Dynamics of 10 random pure initial states 
of qubit Q pictured in the Bloch sphere. Each dot shows the qubit 
state (evolved from the corresponding input state) at a given instant 
of time. Different colors stand for different initial states. Time in- 
creases as the dots move towards point (0, 0, 0), thus showing that Q 
converges towards a state close to a fully mixed one. 



the trace distance V[p^\t)^ p( 2 \t)] at the basis of our cho- 
sen measure of non-Markovianity, one can infer an interesting 
general trend: for values of h/J < 1/2, V generally tends 
to zero at long evolution times. In order to put such an effect 
in context, here we study the long-time distribution of density 
matrices resulting from the evolution of a large sample of ran- 
dom input states according to the channel identified through 
quantum process tomography. Fig. [7] (a) shows the agreement 
between the effective description provided here and the true 
physical effect induced on the state of Q by T. The red dot 
close to the center of the sphere shows the ensemble of 500 
random input states evolved until a final time equal to 100. 
Clearly, there is a single fixed point of the map towards which 



(a) (b) 




FIG. 7: (Color online) (a) Distribution of density matrices evolved 
according to the dynamics induced by a quantum channel corre- 
sponding toh/J= 1/2. We consider a sample of 500 random initial 
states. Regardless of the initial density matrix, the final state is the 
same, thus resulting in a zero-volume state distribution. This is no 
longer the case for h/J = 0.6 [panel (b)], where a finite- volume 
distribution of final states is found. 
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any input state converges after a sufficient amount of time. 

As anticipated, by repeating the above analysis for 
ft/ J > 1/2, we get a distribution of long-time evolved output 
states which is much more widespread, even at very small de- 
viations from the Markovian case. An example is given in 
Fig [11(b) for ft/ J =0.6. At long interaction times, the channel 
does not bring its input states to the same unique state, but a 



finite-dimensional subspace of asymptotic states exists. 

Going back to our attempt of describing the dynamics as 
an effective channel, we have considered the case of a qubit 
evolving under a Markovian channel resulting from the action 
of a generalized amplitude damping channel and an indepen- 
dent dephasing mechanism. The output state for such a pro- 
cess can be written as 



W0)e- 2 ^ +1 ) £ + f^[l - e - 2 7(2A+i)*] 
v pio(0)e-[ 2f +^ +1 )] £ 



An(0)e-[ 2f +^ +1 )]* N 
^[^T-Poo(0)]e- 2 ^ 2 ^ y 



(30) 



r 



with fi being the mean occupation number of the thermal envi- 
ronment responsible for the generalized amplitude damping, 
7 the rate of amplitude damping and F the rate of dephas- 
ing. From this expression, one can determine the matrix of 
the corresponding process, to be compared to the matrix of 
the process associated with the qubit evolution induced by its 
coupling to the chain determined by the quantum process to- 
mography machinery. Calling x(t) suc h matrix evaluated at 
time t and Xc the matrix of the process embodied by Eq. d30b . 
the similarity between such two processes, as time passes, is 
determined by the process fidelity 



Tr(VvWXcVxW) 



(31) 



In calculating this, we have considered scaled parameters 
fi = jit, 7 = jt and F = Ft and optimized F p against them. 
Therefore, time-dependent values of such quantities are re- 
tained at each instance of Jt for the maximization of the pro- 
cess fidelity, which is always at least 90%, as shown in Fig. [8] 
This gives strong numerical indications that our process is an 
homogenizing, time-dependent Markovian channel. 



Markovianity is exactly null, so that the qubit evolution corre- 
sponds to a Markovian dynamical map. By using techniques 
typical of quantum process tomography, we have been able to 
shed light onto the reasons behind the forgetful nature of such 
peculiar working point, showing, in particular, that the zero- 
measure dynamical map has a single fixed point correspond- 
ing to an almost completely mixed qubit state. We have found 
an excellent agreement between the effective Markovian map 
and a generalized amplitude damping channel with an addi- 
tional (independent) dephasing mechanism. The comparison 
between the two descriptions allows us to determine the be- 
havior of the effective dephasing and damping rates, which are 
time-dependent yet positive, in agreement with the Markovian 
nature of the dynamics at hand yfl. 

The richness of the qubit open-system dynamics is very in- 
teresting. We believe that a more extensive exploration of the 
possibilities offered by the tunability of the degree of Marko- 
vianity in such system should be performed so as to under- 
stand if a proper and arbitrary "guidance" of the qubit state 
via the sole manipulation of the properties of the environmen- 
tal chain is in order. This topic is the focus of our currently 
ongoing work which will be reported elsewhere ll2lll . 



III. CONCLUSIONS 



We have studied the dynamics of a qubit interacting 
via energy-exchange mechanisms with a spin chain un- 
der the viewpoint of a recently proposed measure for non- 
Markovianity |3]. We have provided an extensive charac- 
terization of the effective open-system evolution experienced 
by the qubit, often giving fully analytical expressions for the 
measure of non-Markovianity and an operative recipe for its 
experimental determination under certain conditions. Our in- 
vestigation has allowed us to reveal various features of the 
non-Markovianity measure and to relate them to general prop- 
erties of the overall Q — F system. In particular, we have 
focused on the existence of unexpected points in the param- 
eter space of the model, where the chosen measure of non- 
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FIG. 8: (Color online) (a): Optimized quantum process fidelity between a generalized amplitude damping channel and the dynamical map 
here under scrutiny. The process fidelity has been maximized with respect to the parameters of the generalized amplitude damping channel at 
each instant of time. In panel (b) and (c) we show the time-behavior of \i and I\ An analogous behavior is found for 7. The time-dependence 
of the last two parameters, and in particular their positivity, are fully in line with the Markovian nature of the qubit effective dynamics. 
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